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When high-energy cosmic rays (7's, protons, or heavy nuclei) impinge onto the 
Earth's atmosphere, they interact at high altitude with the air nuclei as targets. 
By repeated interaction of the secondaries an 'extensive air shower' (EAS) is gen- 
erated with huge particle numbers in the maximum of the shower development. 
Such cascades are quantitatively simulated by the Monte Carlo computer program 
CORSIKA. The most important uncertainties in simulations arise from modeling 
of high-energy hadronic interactions: a) The inelastic hadron-air cross sections, 
b) The energies occurring in EAS may extend far above the energies available in 
man-made accelerators, and when extrapolating towards higher energies one has 
to rely on theoretical guidelines, c) In collider experiments which are used to ad- 
just the interaction models the very forward particles are not accessible, but just 
those particles carry most of the hadronic energy, and in the EAS development 
they transport a large energy fraction down into the atmosphere. 
CORSIKA is coupled alternatively with 6 high-energy hadronic interaction codes 
(DPMJET, HDPM, neXus, QGSJET, SIBYLL, VENUS). The influence of those 
interaction models on observables of simulated EAS is discussed. 



1 Introduction 

CORSIKA (COsniic Ray Simulation for KAscade) is a detailed Monte Carlo 
program to study the evolution of extensive air showers (EAS) in the atmo- 
sphere initiated by various cosmic ray particles. Originally, it was designed 
to perform simulations for the KASCADE experimenild at Karlsruhe and has 
been refined considerably since its first version in 1989. Meanwhile, it has de- 
veloped into a tool that is used for more than 30 experiments worldwide. The 
prediction of particle energy spectra, densities or arrival times to be observed 
in EAS experiments is a well suited application of CORSIKA. A detailed de- 
scription of the physics incorporated in CORSIKA is given in Ref.H, technical 
details on the handling of the program are described elsewhere.0 

2 EAS Environment Parameters and Particle Transport 

To simulate the evolution of EAS global parameters have to be specified: 
The Earth's magnetic field affecting the movement of charged particles as 
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well as the atmospheric model to be employed in the simulation depend on 
the geographic location. CORSIKA provides several atmospheric parameter 
sets covering the complete climatical and seasonal influence from tropical, 
subtropical, and mid-latitude regions to the South pole (4 seasons). The 
composition of air is adopted to 78.1% N2, 21.0% O2 and 0.9% Ar (volume 
fractions). 

Within CORSIKA, various particles are followed: Besides 7-rays the lep- 
tons e^ and (fc and i^^ optionally), the mesons 7r°, n^, Kg^j^, K='= and 
77, the nucleons, the strange baryons with strangeness \S\ < 3, and the cor- 
responding anti-baryons as well as nuclei with ^ < 60 are treated. Mesonic 
(p, a;,K*) and baryonic (A) resonances are decaying without transport. 

For each particle its transportation range is estimated. For instable par- 
ticles both ranges for interaction (determined by cross section) and decay 
(limited by lifetime) are evaluated independently and the shorter one deter- 
mines the fate of the particle at the end of its range. In decays, all branches 
down to the 1 % level are considered with correct kinematics in the 3-body 
decays. In the range determination of decaying charged particles, the ion- 
ization energy loss is considered, which especially affects muons at energies 
below « 10 GeV because of their long lifetime and low interaction cross sec- 
tion. This treatment slightly favors the decay of charged pions - the main 
source of muons - at the expense of pion-induced interactions, dependent on 
energy and height of pion origin. During transport, the deflection of charged 
particles by the Earth's magnetic field is considered. 

In CORSIKA, j^lectrons and gammas are treated in a tailored version 
of the EGS4 cpdecl and/or, less detailed but much faster, by analytical 
NKG-formulas.El The electromagnetic interactions are believed to be treated 
correctly even at the highest energies by taking into account the Landau- 
Pomeranchuk-Migdal effect. Cherenkov photons may be generated option- 
ally. Hadronic iijteractions with energies E\ah < 80 GeV are modeled by the 
GHEISHA codefl or, alternatively, by the UrQMD model.Ll 

3 Hadronic Interaction Models 

The largest uncertainties in numerical simulation of EAS with primary ener- 
gies above some TeV are induced by the models which describe the hadronic 
interactions. Especially when extrapolating to the highest energies, where the 
collision energies exceed those accessible by accelerators, one has to rely on 
theoretical guidelines to describe the interactions. Additional uncertainties 
stem from the fact that just those interaction products emitted at small an- 
gles into the extreme forward direction carry away the largest energy fraction. 
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Table 1. Essential features of hadronic interaction models. 



Model 


VENUS 


neXus 


QGSJET 


DPMJET 


SIBYLL 


HDPM 


Version 


4.12 


2 




II.4/II.5 


1.6/2.1 




Gribov-Rcggc 


+ 


+ 


+ 


+ 






Mini- Jets 




+ 


+ 


+ 


+ 




Sec. Interactions 


+ 


+ 










N-N interaction 


+ 


+ 


+ 


+ 






Superposition 










+ 


+ 


Max. Energy (GeV) 


2-10^ 


2-10** 


>10" 


>10" 


>10" 


10** 


Memory (Mbyte) 


21 


101 


10 


52 


9 


8 


CPU-Time^ (min) 


4.5 


fSlOO 


1.0 


3.5 


0.75 


1.0 



for showers with primary p, Eo = 10^^ eV, vertical, B^, B^, > 0.3 GeV, 110 m a.s.l., 
NKG option, DEC 3000/600 AXP (175 MHz) 



but in collider experiments those particles disappear in the beam pipe without 
being observable. In the development of EAS, such particles are responsible 
for transporting the energy down into the atmosphere. 

To study the influence of models on the uncertainties of EAS observables 
and on_their cormlations, 6 different hadxpnic interaction jix)des DPMJETB, 
HDPMH, neXusEI, QGSJETy, SIBYLlMia, and VENUsN have been cou- 
pled with CORSIKA. They describe the hadronic interactions at energies 
E\a,h > 80 GeV. Their basic features are summarized in Table ||. 



4 Cross Sections 



The longitudinal development of hadronic EAS depends crucially on the in- 
elastic hadron-air cross section. Lower cross sections elongate, higher ones 
shorten the longitudinal development. Fortunately, the situation has improved 
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in the last 3 years, as demonstrated in Fig. 1. The cross section differences 
between the models have shrunk from 80 mb to today 20 mb in the region 
at a few PeV. Extrapolations to the highest energies above 10^" CleiZ. become 
much more realistic to make predictions for the Auger experiment .E3 

5 Comparison with Collider Data 

All models are adjusted to experimental data wherever available. Especially 
pp collider data should be reproduced with respect to e.g. the increase of 
average charged particle multiplicity with increasing energy, the spread of 
the number of emitted charged particles around the mean value according 
to a negative binomial distribution as shown in Fig. 2a, and the pseudora- 
pidity distribution q£ emitted charged particles shown together with recent 
experimenlaLvalueaHj at i?cm = 630 GeV in Fig. 2b. The more modern 
modelsl3'aE3ll3 displayed here show good agreement with experipacntal data. 
A complete comparison of the older models can be found in RefB3. 

For some parameters, the CORSIKA interaction tests revealed differences 
responsible for deviating properties of FAS simulated with different models. 
The most frequent hadronic interaction within the development of an FAS 
cascade is the collision of a charged pion with a ^^N nucleus. Therefore, 
we compare in Fig. 3a the model predictions of the average charged particle 
multiplicity produced by this reaction as a function of energy. In general, 
the higher the multiplicity, the more energy is dissipated into the production 
of secondaries, which gives rise to a shorter EAS longitudinal distribution. 
The highest multiplicities are observed for QGSJET. In those events with 
many secondary particles only a small energy fraction is left for the relic of 
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Figure 2. Model adjustments to experimental valucs£i|j23 of PP collisions: Charged particle 
multiplicities at = 900 GeV (a), Charged particle pseudorapidity distributions at \/s = 
630 GeV (b). 
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Figure 3. Model predictions for hadron-^^N collisions: Average charged particle multiplic- 
ity as a function of energy for 7r--'^*N collisions (a), energy fraction (Feynman-x) distribution 
of the most energetic baryon from p-^*N collisions at -Eiab = GeV (b). 



the projectile. This is demonstrated in Fig. 3b, which shows the Feynman-x 
distribution of the most energetic baryon from p-^'^N colhsions. In the range 
0.1 < Xp < 0.4, corresponding to a moderate energy transfer to secondary 
particles, QGSJET exhibits the lowest rate. 

Of special interest is the behavior of the Feynman-x distribution of Fig. 3b 
for values Xp > 0.85. This region is governed by diffractive interactions. In 
diffractive collisions the projectile looses only a small energy fraction, the 
largest energy portion is transported deeper into the atmosphere with the 
projectile remainder. It must be pointed out that describing correctly the 
diffractive phenomena is of gteat importance for the hadron rates which are 
observed at KASCADE level.Ell The large differences between the models re- 
flect the lack of experimental data in the forward region. 

6 Results of EAS Simulations 

Only new or recently revised modelJi'B0lli arCpCegarded here. A more de- 
tailed comparison of the models is given in Ref.ESnJ. All shower properties 
which can be measured with detectors placed at observation level more or less 
strongly depend on the 'age' of EAS development. As mentioned above, the 
inelastic hadronic cross sections as well as the 'inelasticity' (the fraction of 
energy carried away by secondary particles) of hadronic interactions are the 
essential quantities which determine the elongation (slow aging) or shortening 
(fast aging) of the EAS development. 

In Fig. 4a the number of electrons (inch positrons) is shown as a function 
of atmospheric depth for vertical showers induced by protons or _Fe-nuclei 
with energies of Eq = lO^"* and 10^''^ eV. The number (averaged over 500 
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Figure 4. Longitudinal shower development for proton and iron induced vertical showers 
at primary energies of Eq = 10^* eV and 10^^ eV (a). Depth of shower maximum Xmax as 
a function of energy as predicted bF- QG i SJET a nd two versions of SIBYLL, together with 
measurement points of experiment ^ ^ ■ P^P^P^r'^ (b). 



simulated showers each) of electrons at sea level (1036 g/cm^) differs between 
the models by max. 14 % (proton) resp. 3 % (Fe) at 10^^ eV, becoming even 
smaller at lower energies. The differences have shrunk by more than a factor 
of 3 with respect to a 1997 comparisor£3 of the older models. 

A quantity measured by several experiments is the depth X^ax (expressed 
in g/cm^) of the maximum number of charged particles within the EAS de- 
velopment. The Xniax- value is sensitive to the primary particle type and may 
be used to determine the (energy dependent) cosmic ray mass composition. 
Fig. 4b givesL-a^aurvey of several experiments,jKhicli use different techniques 
(fluorescenceEj'EZl = filled symbols, CherenkovEScjCJ = open symbols) to de- 
termine Xinax at various energies, extending over 5 decades. 

The model predictions of Fig. 4b are derived from averages over 500 verti- 
cal proton-induced showers (resp. 200 _Fe-induced showers) for each of the 13 
equidistant reference energies. The uncertainties of the mean X^ax are dom- 
inated by shower fluctuations and range from 5 g/cm^ at 10^'* eV to 3 g/cm^ 
at 10^° eV for proton EAS and amount to about the half for Fe. The mean 
A^max-values are approximated by quadratic expressions of the form 
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which follow the simulated mean values within the error bars. 

It must be emphasized that none of the models predicts a linear rela- 
tion as is often assumed in oversimplified arguments. The elongation rate 
per energy decade decreases from c. 70 g/cm^ at 10"^'* eV to c. 50 g/cm^ at 
10^'^ eV for proton-induced showers and shows a similar reduction from c. 80 
to c. 60 g/cm^ for i^e-induccd showers. Remarkable is the good agreement 
between QGSJET and SIBYLL 2.1 up to lO^^ eV. At higher energies, SIBYLL 
predicts a larger Xmax-separation between proton and Fe showers. The new 
SIBYLL version 2.1 generally reveals Xmax- values reduced by « 20 g/cm^ 
relative to the older version 1.6. 



7 Conclusion 



In the last 3 years large progress is attained in EAS simulations. By the 
reevaluation of inelastic proton-air cross sections, a considerable agreement is 
achieved now up to the lO'^ GeV range. A clear trend of convergence between 
different hadronic interaction models is obvious for primary energies up to 
the 10^ GeV range. Despite its age, presently QGSJET still shows the best, 
though not in all respccts Satisfying agreement with a variety of experimental 
resultsMEJ The new SIBYLL versionlij and the more modern ideas realized in 
NEXuga should enable a more consisterLt-and reliable extrapolation especially 
to the highest energies. The announcedE3 version III of DPMJET still has to 
demonstrate its improved quality. 
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